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ABSTRACT 


An analytical study was undertaken to determine the 
two-dimensional wall conduction effects in an internally 
finned, rotating heat pipe. The finite element method was 
employed to generate computer results for a copper condenser 
with triangular fins. Heat transfer rates were shown to be 
apemoximately seventy-five percent creater tham that pre- 
dicted by an earlier, one-dimensional analysis. Heat transfer 
Gaees were found to be insensitive to rotational speed and 
fin half-angle. Due to numerical difficulties within the 
finite element program, no data was obtained for the finned, 


stainless steel condenser. 
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Lo. INTRODUCTION 


Peeeetric ROTATING». NON=CAPILLARY HEAT PIPE 

The rotating, non-capillary heat pipe is a self-contained 
device designed to transfer large quantities of heat in a 
rotating system. The heat pipe contains an evaporator, a 
condenser, and a working fluid as shown in Figure 1. The 
working fluid is vaporized by heat addition at the evaporator. 
The vapor flows toward the condenser end as a result of a 
vapor pressure differential. There, the vapor condenses on 
the inner wall, releasing its latent heat of vaporization. 
The condensate then returns to the evaporator where the cycle 
is repeated. 

The difference between the conventional heat pipe and the 
bemecing, non-capillary heat pipe is the condensate return 
iechanism. A conventional heat pipe employs the capillary 
pumping action of an internal porous wick to transport the 
condensate to the evaporator. The disadvantage of this scheme 
is the resistance to condensate flow through the wick. Addi- 
tionally, the presence of nucleate boiling in the evaporator 
could create bubbles inside the wick structure, reducing the 
effectiveness of the pumping action lS In the rotating 
ieee pipe, rotation about its longitudinal axis establishes 
a centrifugal force field which, when combined with a slight 


I Numbers in brackets indicate references listed in the 
Bibliography. 
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mamrer Of the condenser sect#on, forces the condensate toward 
the evaporator. At high rotation rates, the heat transfer 
rate is significantly increased over the rate achieved in 


conventional heat pipes of similar dimensions. 


B. BACKGROUND 

Previous analyses [2-9] of the rotating heat pipe have 
been conducted to determine its performance characteristics. 
These studies determined, among other things, that under 
Meemal circumstances, the heat transfer rate is limited pri- 
marily by the amount of heat which can be transferred through 
the condenser wall. 

In an effort to improve the condenser heat transfer per- 
formance, Schafer [8] estimated and compared the thermal 
resistances across the condensate film, the wall, and the 
outside convective film using Woodard's earlier data [9]. 

He showed that with the use of outside spray cooling and a 
Bamrecopper wall, the controlling resistance is the inside 
film resistance. Schafer suggested that this resistance can 
be reduced (1) by the proper selection of working fluid, and 
(2) through the use of internal axial fins as sketched in 
Paeure Za. 

Maenining axial fins into the inner surface of the con- 
denser would improve heat transfer performance through two 
Meenanisms: a decrease in condensate film thickness near 
the apex of each fin, and an increase in inner wall surface 


area. However, with the addition of fins, the wall thermal 
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mess tamce would also increase due to the increase in 
effective wall thickness. Of course, the use of copper as 
the wall material, because of its high thermal conductivity, 
would minimize the added wall resistance. 

Schafer then developed an analytical model of the intern- 
[my tinned condenser for the case of fins of triangular 
Protizle (Figures 2a, 2b). In the analysis, he assumed one- 
dimensional heat conduction through the wall and assumed a 
linearly decreasing temperature profile from the fin apex 
(point A) to the midpoint of the trough (point C). Results 
Maomethne analysis indicated significantly improved heat 
transfer rates for the finned condenser compared to the smooth 
condenser. He pointed out, however, that two-dimensional 
conduction effects may be significant so that a more detailed 
analysis of the heat conduction problem within the wall was 


necessary. 


Gao Se ECTIVES 

The objectives of this study were, therefore, to (1) 
develop an analytical solution to account for two-dimensional 
conduction effects within the wall for the triangular fin 
meoerle, and (2) compare this solution with the solution 


@enived by Schafer. 
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Pi. IJHEOREILICAL PROGRAM 


A. REVIEW OF SCHAFER ANALYSIS [8] 

Schafer's study considered an internal fin geometry 
represented by Figure 2a. Since the sum of wall thickness 
and fin height was assumed to be much smaller than condenser 
radius, an average radius from any point on the axis to the 
condenser wall was taken to be the distance from the axis to 
a position midway between the base of the fin and the apex 
otetne fin, ri: The centrifugal acceleration term was then 
evaluated at this average radius, Tw’. 

The thin wall assumption also allowed the curvature of 
the wall to be neglected. The problem could then be analyzed 
using a planar conduction heat transfer model as shown in 
Fosemre 2b. Furthermore, because of repetitive geometry the 
problem could be reduced to the simpler geometry shown in 
Figure 3. Cartesian coordinates were then selected, as shown, 
to permit a simplified fin condensation analysis. 

Assumptions in the analysis were based upon the earlier 
work of Ballback [2]. These assumptions were similar to those 
used in the Nusselt analysis of film condensation on a wall. 
The more important of these assumptions are: 

(1) steady state operation, 
(2) film condensation, as opposed to dropwise condensation, 


(3) laminar condensate flow along both the fin and the 
trough, 


(4) static balance of forces within the condensate, 


is 
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(5) one-dimensional conduction heat transfer through the 
film thickness (no convective heat transfer), 


(6) no liquid-vapor interfacial shear stresses, and 
(7) no condensate subcooling. 
Additional assumptions in the analysis are noted later. 

The complete development of Schafer's finned condenser 
heat transfer theory is presented in his thesis [8]. However, 
Since the fluid flow and heat transfer theory was applicable 
to the two-dimensional conduction heat transfer analysis, a 
summary of his work is presented here, prior to the develop- 
ment of the two-dimensional conduction model. 

1. Fin Condensate Mass Flow 

A Nusselt analysis was performed for film condensation 
on the wall of the fin. In this analysis, the fluid momentum 
equation was constructed to account for the effects of condenser 
cone half-angle, 9, fin half-angle, a, and the centrifugal 
pomee field. The pressure differential term was megdected. 
Condensate flow on the fin was assumed to be in the z-direction 
OMey . 

An expression for fluid velocity, w(y,z), was then 
derived by substituting the shear stress equation for laminar 
flow, ce > Ve ~ » into the momentum equation and solving 
for w(y,z). Average velocity, w(z), was then found by inte- 
Cmeeing w(y,z) over y from y=0 to y=46(z), the local film 
thickness. This average velocity was then placed in the fin 
mass flow equation to yield an expression for fluid mass flow 
rate along the fin per unit condenser length as a function 
Oise Z), 
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a * 





pw (R, +x sing) cos¢cos a 6(z) 


—— (1) 


Ds roush Condensate Mass Flow Equation 


Fluid flow within the trough was assumed to be one- 
dimensional in the x-direction and of constant velocity across 
the width of the trough. Therefore, condensate flowing from 
the fins into the trough was assumed to accelerate immediately 
to the velocity of the trough fluid. Cross-sectional area of 
tie trough condensate region was twice the area shown in 
Figure 4. 

The x-direction momentum equation for an infinitesimal 
fluid element in the trough was constructed with the pressure 
(feeeerential term neglected. The trough mass flow expression 
Woemernen developed in a manner similar to the development of 
the fin mass flow equation (1). The average trough velocity 
Was derived and substituted into the trough mass flow equation 


to yield total mass flow through the trough region, 


ppw(R +x sin 6) 8* (x) Sin 


eee) aa su, 


(een) + Sin a cose et. (Z;) 


where the expression in the second parenthesis is the cross- 
sectional area of the trough condensate region. 

Total condensate mass flow was identical to trough 
condensate mass flow since all condensate eventually flowed 


to the trough. Therefore, once the total mass flow rate at 
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any position x was known, trough condensate thickness, 6*(x), 
could be determined using equation (2). 
3. Condenser Heat Transfer Determination 

Schafer assumed zero heat flux boundary conditions on 
both sides of the wall section shown in Figure 4. Temperature 
at the fin apex was taken to be fluid saturation temperature 
since film thickness at this position was assumed to be zero. 

few etc Mmpescea « linearly deercasing temperauure on 
eiceanner wall surface fromthe apex of the fin to tive mid- 


channel position. The surface temperature expression became 
Ts - T,@4) = 82, (3) 


where 8 was an unknown constant. This choice of temperature 
distribution was arbitrary but allowed a simplified one- 
dimensional conduction analysis of the condenser wall. He 
then imposed the restriction that no heat transfer would 
occur through section BE of the wall (Figure 3). All heat 
transfer through the fin surface, thus, flowed out of the 
wall through section FE. 

Heat transfer through the fin was equated with one- 
dimensional heat transfer through a rectangular geometry of 
the same width as section FE, and of height, t*, equal to t 
plus one-half the fin height (Figure 4). Temperature of the 
upper surface of the rectangle was set equal to the tempera- 
ture at a point midway between positions A and B. The 
expression for heat transfer rate per unit condenser length 


through the wall and outer convective thermal resistance 
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of the fin became 


»  ((T, - 2) - TL] e, 


1 tc 1 ae (4) 
CaS ae) 


+ 
W out 





It was then desired to develop an expression for 
Q (2) in terms of the heat transfer rate through the fin 
condensate film. First, the equation for fin condensate 
film thickness was derived by differentiating (1) with respect 
to z and substituting the result into the differentiated 


equation for mass flow rate per unit condenser length, 


5 r dQ. (z fel We S32) 

aM, (z) = SAte) 2 2 fx, (Sul) az. (5) 
fg fg 

The expression produced was algebraically rearranged to form 


a differential equation in 6(2Z), 


kei Tra) dz 
ez)a6(z) = ie aa aa (6) 


2 2 : 
Pe w (R, + X Sanco Ngo cos ¢ cosa 


Equation (3) was then substituted in (6) and the result 
iiieeerated over z from z = 0 to z to give an equation for 
6(z), 
| Zul B Ue 1 
i. (eee oe, 
og w'(R, +X sin >)he, cos ¢ cosa 
The value for fin heat transfer rate per unit condenser 


length was then evaluated using 
* = ’ * 
Q, (z*) = hg, M (z*) (8) 


Jey 





Wieme 2° is defined in Figure 4. This expression was set 


equal to (4) to yield a transcendental equation in B, 





[(T, - oe - halle Gs 
@. oS) 
ee | Sie 


V/; 
(2K -8)*/ (9 -w) 2 [(R,+xX.sin 6) he, cos $ cos a] "aah 


, 
Sug (9) 
where E, was the width of the base of the fin. The expressions 
on the left and right hand sides of this equation were equiva- 
lent only when all boundary conditions were satisfied. The 
task was then to solve iteratively for 8 until equation (9) 
was satisfied. Q, (2*) was then determined in equation (4) 
using the final iterated value of 8. 


Heat transfer rate from the trough region was evaluated 


us ing 
: (T. - T ) 
Omix)ye=n ks 2 ee >; (10) 
6 * (x) 
where Tove was Calculated as an average temperature of the 


Geawen wall and € was the total trough Width (Figure 4). 
Peetvowenarewacscription of the computer program used 

ieecie Schafer analysis is given in Figure 5. The axial 
length of the condenser was first divided into increments 

of equal length, Ax. Solution of the heat transfer problem 
was then started by assuming an initial trough condensate 
thickness, oan = 0.00006752 feet. This initial value was 
based upon a mass flow analysis by Sparrow and Gregg of 


Gondensacwon on a rotating disk [10]. It was assumed in 
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Schafer's analysis that the effect on the Sparrow and Gregg 
solution of enclosing the disk within the condenser would 
be minimal. A study of the effects of a ten-fold increase 
or decrease from the above value of 6*, indicated that the 
initial 6*, assumption would not affect the heat transfer 
solution by more than ten percent. Therefore, the deter- 
mination of an exact value of initial trough condensate 
thiekness was not crucial to the problem. 

Calculation of condensate mass flow contribution from 
the endplate then followed. Then, for each increment Ax, an 


Bz 
fimenmal value of 8B was Calculated by removing the —,2 term 


2 
from the left side of equation (9) and solving directly for 
8. The full equation was then employed to solve iteratively 
for 8. Calculation of incremental fin heat transfer rate 
then followed by multiplying the value obtained in (4) by 
the incremental length, Ax. Incremental trough heat transfer 
rate was next evaluated by multiplying equation (10) by Ax. 
meememental trough and fin heat transfer rates were then 
Summed to yield the total incremental heat transfer rate for 
the geometry of Figure 4. 

Total incremental mass flow rate was evaluated by 
dividing twice the sum of the incremental trough and fin heat 
transfer rates by the latent heat of vaporization of the fluid. 
As the incrementation progressed, incremental mass flow rates 
were added to mass flow from preceding increments to yield 
M (x). Equation (2) was then solved for 6*(x) with the 


EOL 


aidsor a polynomial rootfinder subroutine. The new value of 
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6*(x) was then used to begin the solution for the next 
increment. When the total length of the condenser had been 


traversed, total heat transfer rate was evaluated. 


: 100 : : 
Qtot 7 Nein ra Ace (2 eet OnGa)eI - Gli} 


B. TWO-DIMENSIONAL HEAT TRANSFER SOLUTION 
mee ncroduc tion 

In this study, two-dimensional conduction through the 
condenser wall was considered in the heat transfer analysis 
using a finite element solution. The finite element method 
was chosen over the method of finite differences because of 
considerations involving wall geometry and computer utiliza- 
tion time. Myers [11] has shown that when these methods are 
compared for the case of two-dimensional conduction in a 
rectangular block with heat generation, the finite element 
method more closely approximates the exact solution when the 
number of nodal points considered in each method are equal. 
This is particularly true when the nodal point meshes are 
coarse. 

To achieve the same accuracy as with the finite element 
method, the finite difference method must use a finer nodal 
mesh, especially for non-rectangular geometries. Additionally, 
when fins of different shape are considered, the mesh of the 
momrce difference model must be completely reconstructed for 
each new shape. In the finite element program, however, a 


Simple change in nodal coordinates is all that is necessary 
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to accommodate the change in fin shape. Since future studies 
may analyze the use of. varied fin shapes with the rotating 
heat pipe, it was desired to incorporate the more flexible 
finite element method in the computer program to avoid having 
to rearrange the program for each fin shape. 

While the accuracy of the finite element method is 
Siieerior, the finite difference method requires less computer 
ierbezation tame for the same number of nodal points. If a 
finite difference model of finer mesh was chosen, the decrease 
mmconouter time due to the method might be offset by an -in- 
crease in time required for the increased number of nodes. 
Thesoverall improvement in computer utilization time might, 
therefore, be minimal. As a result, in this study, the choice 
of method was not based upon computer time requirements. Also, 
a finite element computer program [12] was available which 
could be adapted for the purposes of this analysis. For the 
reasons outlined above, therefore, the finite element method 
was selected. 

The steady state differential equation for conduction 


heat transfer is given in Arpaci's text [13] as 


oT e) w 
apy oy ee <a (12) 


Boundary conditions are normally of two kinds; those of 


specified temperature on a boundary, 
T = Ty , (boundary b only) (ars) 


and those of specified heat flux or convective condition 
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om the boundary, 
dT 
Ky aly dy (ky az 4 2p 


+ {h(T - ere + oo 0, (boundary b only) 
(de 
where cy and ¥ are direction cosines of the outward 
normal to the boundary's surface. 
fee ine Fingwee Element Method 
Many references including [11], [14], and [15] develop 
the concept of using material elements of finite dimensions 
Within a structure to model the differential equation and its 
boundary conditions. The steady-state field equation, (12), 
may be discretized by considering that the temperature at any 
Pomitewithin an element is a function of the temperatures of 


specified nodes on that element, 


CNT =e (15) 


where <N> is the row matrix shape function of position and 
{T}° is the nodal temperature vector for the element. The 
field equation and boundary equations for each element may 

be mathematically rearranged and integrated to yield a matrix 


Syoeem OL equations for that element, 


e 


Gl) lS Sets ater (16) 


where 


[H] © 


it) 


f CE) + # en >) )dz dy , ee) 


Se 


and 


ae 











or oe - (fq {N} de), 


Ce 


+ (Ug, (Nb KN) de), - ALT}, op. (18) 


Cc 


Equation (18) is applied only over relevant boundaries. The 
element matrix equations are then assembled into a master 
meerix equation for the entire structure. Solution by 
numerical methods then follows. 

The comvective boumdary conditions were applved in 
tert inite element solution of the finned condenser using a 
lumped convective condition formulation. Average nodal con- 


vective coefficients, h were determined for each node 


avg’ 
subject to a convective environment. The nodal have values 
were then considered to be constant over the boundary region 
associated with the node. The method of determination of 
these coefficients is described in a later section. 

Equation (18) was, thus, reduced to the following 


form, 


ce ec | .e 
(Fly = (hay, Th - (hay, Tah Jp. (19) 


This lumping procedure reduced the complexity of the convective 
condition calculations and did not seriously affect the values 
Within the master matrix, [H]. 

Reference [15] presents a development of the 
meparametric finite element concept. This concept allows 
Beieasolution of eqtt@taons (12), (15), and (14) for elements 
of non-uniform shape. The application of the method of finite 
elements is, therefore, extended to a wider range of physical 
problems. 
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The eight-nodal finite element mesh used in the 
two-dimensional analysis of the axially finned condenser is 
shown in Figure 6. A computer program for the case of two- 
dimensional conduction was developed using a modified version 
of the three-dimensional, isoparametric finite element program 
assembled by Lew [12]. 

oer waest Cases 

Betoremthe. finate .lemnent.soimeron, was,applaed, to. the 
condenser problem, a series of test cases were analyzed. 
These cases allowed a comparison of analytical solutions to 
results obtained using the finite element program. The two- 
dimensional heat conduction problems studied were 

eee rectangular block with specified surface temperature, 


(2) a rectangular block in a uniformly convective environ- 
ment, 


Meee rectangular block with a triangular fin in a uniformly 
comvective environment, 


(4) a rectangular block with a convex parabolic fin ina 
uniformly convective environment, and 


(5) a rectangular block with a concave parabolic fin in a 
uniformly convective environment. 


Results of the test cases are presented in Figures 7-11 with 
values obtained using the finite element program placed in 
brackets. For the first two cases, the heat transfer rates 
Calculated with the aid of the finite element computer program 
agreed well with the exact, two-dimensional solutions. 
Although finite element heat transfer solutions were rela- 
tively accurate, nodal temperature values were slightly 


different from the exact values. This temperature discrepancy 
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led to numerical difficulties when the finite element solution 
Mas applied to the stainless steel condenser. These diffi- 
culties will be discussed more fully in Chapter III. 

For the finned block cases, analytical expressions 
were derived for heat transfer rates through extended surfaces 
using the one-dimensional fin analysis as described in [16]. 
Heat transfer rates calculated by the finite element program 
for this uniformly convective environment did not vary sig- 
meerecancly from the one-dimensional results. Samce test case 
results from the finite element program agreed favorably with 
the derived results, the program was assumed to be operating 
properly and was applied to the condenser heat transfer 
problem. 

4. Condenser Boundary Condition Application 

The axially finned condenser problem was analyzed 
iyeeceamsidering the geometry of Figure 6. Boundary conditions 
mer the problem were considered to be the same as those in 
the Schafer analysis with the exception of the inner surface 
temperature distribution. Instead of using a linearly decreas- 
ing temperature profile along the wall from the apex of the 
fm to the midpoint of the channel, a parabolic distribwtion 
was assumed and used to calculate average nodal convective 
@eerftircients. 

The Nusselt analysis of the fin condensation problem 
yielded a differential equation in 6(z), fin condensate 
thickness, which is identical to equation (6) used in Schafer's 


analysis. 
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Consistent with the eight-nodal finite element 
concept, a parabolic temperature distribution was then 


assumed to exist along the fin surface using the boundary 


condition, T = ie atez =—0, simce G00) = 0. Thepmerore, 
az) = az" 1) anaes Wee (20) 
Ong 
Tas Eee) — = Da. (21) 


where a, and D, were unknown constants. 
When (21) was substituted into equation (6) and the 
result integrated from z = 0 to z, the expression for 6(z) 
became 
3 2 L/ 
Sion ( - ae = bee) u ‘ 
é(z) = (aE ee (am 


Og w(R + x sin 6) hee cos } cosa 


The local convective heat transfer coefficient, assuming 
one-dimensional heat conduction through the thin condensate 
film, was then 

‘1s 


Kee To eee Said) he, cosd cos a 


Pein 3 i ff x z2 
due (- a - oF a ) 

(23) 
and was applicable from z = 0 to z = z* = a o°7 ees OC 
(Figure 6). 

Across the width of the trough, from z = Zz, to the 
mid-channel position, 


Ke 
n(X) trough 7 Seco) * a 
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6*(x) was considered constant across the trough width. Based 
meen a check of the»ewall temperatures generated from Schafer's 
computer program, the fluid temperature difference across the 
width of the trough was predicted to be small. Fluid thermal 
conductivity values were, therefore, assumed to remain 


relatively constant. As a result, h(x) was assumed to 


trough 
be constant across the trough width and was calculated using 
equation (24) and a representative average inner wall temper- 
ature to determine Ke. 

An h(z) temperature distribution was then needed 
between z = z* and Zz = Zoos Since the thickness of condensate 
film increased linearly from the value at z* to the value at 
Zo the convective coefficient was assumed to decrease linearly 
from h(z*) ¢.. to h(X) enough’ 

Other boundary conditions were then applied in 
accordance with the assumptions of Schafer, namely, zero heat 
flux through the sides of the geometry, T = Te at the fin 


apex, and constant convective coefficient, h along the 


out’ 
outside surface of the wall. Boundary conditions for the 
mee iem were, thus, completely specified. 

In his expression for trough condensate cross-sectional 
area, Schafer assumed that the surfaces of the condensate 
mhews from the fin and the trough were coincident at a position 
directly above the trough-fin wall juncture (Figure 4). For 
small fin film thicknesses and other than small fin angles, 


this assumption was not accurate. In this study, it was 


assumed that fin condensate thickness was much smaller than 


cy 
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trough condensate thickness. The trough condensate region 
could, therefore, be more accurately described as shown in 


Figure 6. The modified trough mass flow equation became 


Pe Go GR, ex sin 9)6*(x)}? sin $6 
meet «Ge + 6 A(x) a). 
(25) 
It should be noted that this equation differs from equation 
(2) only in the expression for cross-sectional area of the 
Peete condensate region. The change in cross-sectional areas 
involved the term, 1/cos* 4, which for large fin angles, could 
make a significant difference in cross-sectional areas of the 
two equations. 
5. Condenser Heat Transfer Solution 
Within the computer program, the heat transfer rate 
for the entire condenser was determined by dividing its axial 
length into one hundred equal length increments of length, 4x, 
Soevaine for the incremental heat transfer rate using an iter- 
meive procedure, and summing those incremental rates to yield 
a total heat transfer rate. A flowchart description of the 
eemoucer program used in this analysis 15 given in Figure 12. 
At each increment, constants a, and b, were determined 
using a Lagrangian interpolation fit of the temperatures 
calculated from the previous iterate at nodes 1, 3, and 6. 
hem the first increment, a and by were initially calculated 
from an initial guess of the temperatures at those nodes. 
Convective coefficients at several locations along the fin 


were determined using equation (23) and the linearly decreasing 
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imezyeassumption from z = z* to z = ie Next, average nodal 


values, hy avg’ Ns avg? and the average h of the left hand 
boundary region associated with node 6, (Figure 13), were then 
calculated using the position dependent values of h(Z)es, in 


a Simpson's Rule approximation. h(x) was then calculated 


trough 
from equation (24) using the §* value calculated from the 
previous increment. This value was used as the average 
Convective coefficient for nodes 9 and 12 and the right hand 
boundary region of node 6. The average convective coefficients 
Sreewracted for the left and right hand boundary regions of 

node 6 were then weighted by their respective region lengths. 
The two values were then summed to yield the average convective 


coefficient for node 6, h The remaining boundary condi- 


Savg- 
tions were then applied and incremental heat transfer rate 
femeuieee Condenser length, - and nodal temperature values 
Were generated from the finite element subroutines. 

If the resultant heat transfer rate converged to 


Within the criteria of 0.05 percent change from the value of 


Ene previous iterate, 


Q, - Q 

moo, (26) 
Or 
j 


the incremental problem was deemed solved and Q; set equal to 


the final iterate, Je 


criteria was not met, the newly calculated nodal temperature 


If, on the other hand, the convergence 


values were used to establish a new set of inner wall 


convective coefficients. After the remaining boundary 
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conditions had been applied, the finite element solution was 
repeated until heat transfer rate converged to within the 
Seeempeaple 0.05 percent tolerance. Because Q, was initial- 
ezemewet zero for each increment, the itemation procedure was 
mepeated at least twice for every increment. 


Incremental mass flow rate was then calculated using 


: GF 
AM = 2 —+ Ax. (C20 
tot, Neg 


This value was then added to the mass flow rate from previous 
increments. The subsequent interval's trough condensate 
thickness was then calculated using a polynomial root-finder 
subroutine to solve equation (25). The incrementation 
continued until the entire length of the condenser had been 
traversed. Incremental heat transfer rates were than summed 


aiaearctotal heat transfer rate calculated from (3). 


‘ 100 =: 
Qrot ~ 2 ae & (Q; a), (28) 
ibalb 
where Nein was the number of fins along the inner wall 
circumference. 


To start the incrementation, two initial assumptions 
were made. The initial increment's value of 6* = Se was 
meen trom the analysis of Schafer as 0.00006752 feet. For 
reasons explained in the previous sub-chapter, a large variance 


Maeche accuracy of this value would not significantly affect 


the heat transfer results obtained from the computer progran. 





Additionally, at the first increment, initial guesses 
of fin temperatures at nodes 1, 3, and 6 were made in order 
to begin the solution. As an example, the temperatures used 
Goebpecum the solution for the copper condenser with fluid 
Saturation temperature of 200°F were i = 200°F, ip = 195°F, 
and A SeeliOd° F. 

For the coppes: cenmdenser, thes use of widely variant 
initial nodal temperature values, 1S 2002, iw = 191°F, 
and = 182°F had little effect on convergence to the heat 
transfer solution for the first increment. However, as will 
be discussed in Chapter III, the initial temperature guess, 
coupled with the slight inaccuracy of the finite element 
nodal temperature solution produced numerical difficulties 


within the computer program for the stainless steel condenser. 





Pie) Petiiso LON, OF RESULTS 


Computer results for the finned copper condenser were 
@enerated for three outside convective coefficients, three 
fin angles, and five rotation rates. The working fluid used 
throughout the analysis was water. Representative results 
are presented in Figures 14 and 16, 17, 18. 

Figure 14 compares the heat transfer results of this 
Study with the results of the analysis of Schafer for given 
Operating conditions. It can be seen that the two-dimensional 
Gonductwon analysis allowing a parabolic temperature distri- 
bution along the fin predicted heat transfer rates for the 
finned condenser which are approximately seventy-five percent 
greater than the rates predicted with the one-dimensional 
model with linear temperature profile along the wall. This 
rather significant change was unexpected, and may be some- 
what in error due to numerical difficulties in the method. 
These numerical difficulties arose due to the boundary 
condition at node 1. The finite element boundary conditions 
ame established for each node such that the condition is 
assumed to exist over the entire boundary region of the node. 
Therefore, the zero flux boundary condition could not be 
applied over the whole length of the left and right boundaries 
of the condenser geometry. This problem existed for nodes l, 
@eell, sand 12, but owas most critical for node 1 since. it was 


estimated that twenty to thirty percent of the total heat 
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flux through the inner wall passed through the region 
associated with this node. 

An error of fifty percent in the calculation of node 1 
heat transfer rate might then give a total heat transfer 
error of as high as fifteen percent. Two possible solutions 
to the boundary condition problem merit attention. The first 
is to increase the number of elements used in the solution of 
the geometry of Figure 6. However, a boundary region associ- 
ated with node 1 would still exist on the left hand side of 
the geometry although reduced in size. The effects of the 
boundary condition problem would, therefore, be reduced but 
not eliminated. 

The second solution would be to consider a wall geometry 
as shown in Figure 15. The zero flux boundary condition 
would then be forced along the boundary in question due to 
the symmetry of the problem. The use of one of these alter- 
nate solution methods might contain the added benefit of 
eliminating the numerical difficulties associated with the 
Stainless steel condenser solution. 

The heat transfer rate solution is generally affected in 
the same manner as with the analysis of Schafer when Bess 
Ww, and a were varied. In Figure 16 the influence of outside 


convective coefficient, h 1s shown. Figure 17 illustrates 


eve ¢ 
how the rate of heat transfer was relatively unaffected as 
fin half angle, a, was varied. In addition, heat transfer 


mate was found not to vary significantly with rotation rate 


as shown in Figure 18. This result concurred with the results 


o5 


obtained earlier by Schafer [8]. Schafer suggested that the 
memvariance in heat transfer rate with the rotation rate was 
due to the cancelling effects of decreasing film thickness 
and decreasing temperature difference across the film 
Babekness as rotation rate is increased. | 

Because of numerical difficulties in the program, no 
mesmlts were generated for the staambess steel condenser. 
With the solution of the first and second test cases, the 
nodal temperatures calculated were accurate only to within 
Beomor three degrees. In the condenser problem, this slight 
inaccuracy resulted in nodal temperatures along the fin wall 
which were slightly above the fluid saturation temperature. 
Computer output indicated that the wall temperature at node 3 
was approximately 0.1°F above the saturation temperature when 
ta 200°F. As a result, the quantity within the brackets of 
equation (23) was a negative value. Error messages were 
generated when the computer attempted to evaluate the fourth 


root of this negative number. 
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iv. GONCLUSTONS 


lee the addition of internal axial fins to the condenser 
surface of a rotating heat pipe significantly improves heat 
@eamster performance. 

2. Analysis of the condenser heat transfer problem using 
a two-dimensional wall conduction model and a parabolic 
@emperature distribution along the fin yielded heat transfer 
Hates about seventy-five percent above those predicted by the 
one-dimensional, linear temperature distribution analysis of 
schafer [8]. 

3. For the geometrical cases studied in this thesis, the 
heat transfer rates calculated from the two-dimensional 
analysis are insensitive to changes in rotational speed and 
fin half-angle. 

fee boundary Condition problems produced with the finite 
emement Solution for the copper condenser were the probable 
cause of solution inaccuracy on the order of fifteen percent. 

5. Slight inaccuracy from the finite element calculation 
Of nodal temperature values led to numerical difficulties in 


the heat transfer solution of the stainless steel condenser. 
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V. RECOMMENDATIONS FOR FURTHER STUDY 


It is suggested that the following topics be considered 
for future study of the internally finned, rotating heat 
Pplpe: 

l. Re-examination of the finite element analysis using 
Siemer an increased number of elements for the geometry 
considered in this study or the mirror image geometry pre- 
sented in Figure 15S. 

2. Use of twelve-nodal elements in the finite element 
a@alysis to allow the assumption of a cubic inner wall temper- 
femee protile in order towmore closely approximate the actual 
temperature distribution. 

Seeeenalysis of film condensation on a variety of fin 
Shapes, and application of the solution to the two-dimensional 
model to determine the shape permitting optimum heat transfer 
hate’. 

een experimental program in order to compare theoretical 
and experimental heat transfer rates of the internally finned 


condenser. 
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Evaporator Condenser 


Figure 1. Rotating Heat Pipe Components 





Figure 2a. Internally Finned Condenser Geometry 





Peeure 2b. Panar Comdwetion Simplification 
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Figure 4. Condenser Geometry Considered in Analysis 
of senmarcr [hs] 
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Figure 7. Test Case (1) 
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Figure 10. Test Case (4) 
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Saturation Temperature 
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Figure 15. Alternate Condenser Wall Geometry 
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